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Abstract: We study the effect of a veto on additional jets in the rapidity region 

r^i between a pair of high transverse momentum jets at the LHC. We aim to sum the 

most important logarithms in the ratio of the jet transverse momentum to the veto 

scale and to that end we attempt to assess the significance of the super-leading log- 

^ arithms that appear at high orders in the perturbative expansion. We also compare 

our results to those of H.ERWIG+ + , in an attempt to ascertain the accuracy of the 

angular ordered parton shower. We find that there are large corrections that arise for 

O large enough jet transverse momenta as a consequence of Coulomb gluon exchanges. 
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1. Introduction 

Dijet production, with a veto on the emission of additional radiation in the inter-jet 
region has been widely studied in electron-proton collisions at HERA [1-3] and in pp 
collisions at the Tevatron [4,5], and theoretical calculations have been compared to 
data [6-10]. We shall refer generically to the "gaps between jets" process, although 
it is to be remembered that the veto scale can certainly be large and so a "gap" is 
simply a region of limited hadronic activity. Indeed we shall require the veto scale 
to be large in this paper in order to permit a calculation using perturbative QCD. 

Gaps between jets is an interesting process to study at the LHC. It is also a pure 
QCD process, hence the cross-section is large and studies can be performed with early 
data. It is interesting because it allows one to investigate a remarkably diverse range 
of QCD phenomena. These are summarised in Fig. 1, which maps out schematically 
the (L, Y) plane, where L = ln(Q/Q ) (Q is the dijet pr and Qq the veto scale) 
and Y is rapidity separation of the dijets. For instance, the limit of large rapidity 
separation corresponds to the limit of high partonic centre of mass energy and BFKL 
effects are expected to become important [11]. On the other hand one can study the 
limit of emptier gaps, becoming more sensitive to wide-angle soft gluon radiation. 
Furthermore, if one wants to investigate both of these limits simultaneously, then 
the non-forward BFKL equation enters the game [12]. In this paper we concentrate 
on wide-angle soft emissions. 
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Figure 1: A map of the (L, Y) plane. 



Given a hard scattering process, we can study how it is modified by the addition 
of soft radiation. If the observable is inclusive enough, then we have no effects be- 
cause soft contributions cancel when real and virtual corrections are added together, 
as a result of the Bloch-Nordsieck theorem. However, if we restrict the real radiation 
to a corner of the phase space, as happens for the gap cross-section, we encounter 
a miscancellation and we are left with a logarithm of the ratio of the hard scale 
and veto scale. The resummation of wide-angle soft radiation in the gaps between 
jets process was originally performed assuming that the real-virtual cancellation is 
perfect outside the gap, so that one need only consider virtual gluon corrections 
integrated over momenta for which real emissions are forbidden, i.e. over the "in 
gap" region of rapidity and with k T above the veto scale Q [13-15]. However, it 
was later realised [16] that one must also consider secondary emissions, i.e. real 
gluons emitted outside of the gap may emit back into the gap, and these emissions 
too are vetoed above Qo- These give rise to a new tower of logarithms, formally 
as important as the primary emission corrections, known now as non-global loga- 
rithms. In order to obtain the complete tower of leading logarithms it is therefore 
not sufficient to consider only the virtual corrections to the 2 — > 2 scattering process. 
To capture the non-global logarithms, one must also consider all 2 —>■ n processes, 
i.e. n — 2 out-of-gap gluons, dressed with "in-gap" virtual corrections. The colour 
structure quickly becomes intractable and, to date, corrections have been calculated 
numerically to all orders only in the large N c limit [16, 17]. In [18] it was shown that 
the resummation of non-global logarithms in the large iV c limit can be obtained by 
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solving a non-linear evolution equation, which has the same form as the BK equation 
in small- a; physics [19,20] (see also [21]). This remarkable observation has stimulated 
further investigation into a correspondence between soft gluon emission and high 
energy (Regge) scattering [22-24]. 

A different approach was taken in [25] , where the specific case of only one gluon 
emitted outside the gap, dressed to all orders with virtual gluons but keeping the full 
N c structure, was considered. The surprising result of that calculation was the pos- 
sible discovery of a new class of logarithms. These logarithms appear first at order 
a A s relative to the hard process and are formally more important than the primary 
emission logarithms at this order. These "super-leading" logarithms (SLL) are how- 
ever sub-leading in colour, which is perhaps the reason why they were not identified 
in previous studies. Their origin can be traced to a failure of the DGLAP "plus- 
prescription", when the out-of-gap gluon becomes collinear to one of the incoming 
partons. Real and virtual contributions do not cancel as one would expect and one 
is left with an extra logarithm. This miscancellation is caused by the imaginary part 
of loop integrals, induced by Coulomb gluons. More recently the calculation has 
been repeated in a colour basis independent formalism [26] . Such an approach allows 
one to consider, in principle, any number of gluons outside the gap. However, in 
order to perform an actual calculation, one must fix a colour basis and its increasing 
dimensionality hampers the computation. In this paper we limit ourselves to the 
resummation of SLL contributions coming only from the emission of only one gluon 
outside the gap. 

Accurate studies of these effects are important also in relation to other processes, 
in particular the production of a Higgs boson in association with two jets. It is 
well known that this process can occur via gluon-gluon fusion and weak-boson fusion 
(WBF). QCD radiation in the inter-jet region is clearly different in the two cases and, 
in order to enhance the WBF channel, one can put a cut on emission between the 
jets [28,29]. This situation is very closely related to gaps between jets since the Higgs 
carries no colour charge. To date, the impact of a central jet veto has been explored 
in fixed-order calculations [30-32], using parton shower Monte Carlos [33,34], via 
soft gluon resummation [35] and in the framework of high energy factorisation [36] . 
Control over the effects of QCD radiation in Hjj is mandatory if we are accurately 
to extract the coupling of the Higgs to the weak bosons in this channel [37-40]. 

The remainder of this paper is organised as follows. In Section 2 we define the 
observable and the kinematics. In Section 3 we reconsider the known case and sum up 
the subset of leading logarithms that arise as a result of imposing a veto on primary 
emissions into the gap. We refer to this as the "zero gluons outside the gap" cross- 
section and we study the effect of the resummation on the hadronic cross-section at 
the LHC. We compare our result to HerwigH — h and are able to assess the role of 
the Coulomb gluon corrections, which are not present in the angular ordered parton 
shower approach. In Section 4 we perform the resummation of the SLL contributions 



- 3 - 



that come from matrix elements in which there is one gluon emitted outside of the 
gap and we speculate on the size of the missing contributions from the emission of 
greater than one out-of-gap gluon. We finish with a look to the future. 



2. The gaps-between-jets cross-section 

We are interested in dijet production in hadron-hadron collisions: 

hi + h 2 ^ j j + X , 

where we veto the emission of a third jet with transverse momentum bigger than 
Qo in the rapidity region between the two jets. At leading order in perturbation 
theory we have to consider all the possible 2 — > 2 parton scattering processes. The 
kinematics of such sub-processes can be expressed as follows: 



Pa = X a — (1,0,0,1) 

Pb = x 6 — (1,0,0,-1) 

Pc = Q(cosh y c , cos ip, sin <p>, sinh y c ) 

Pd = Q (cosh y d , - cos <p, - sin ip, sinh y d ) , (2.1) 

where 

x a ,b = ^(e^ + e^) , p=^. (2.2) 



The triple differential cross-section is 



d 3 <7 g 4 J2 fa( X a,Q 2 )fb(x b ,Q 2 ) 



dy c dy d dQ 2 lQnS 2 ^ 

a,o,c,a 

|2 1 1 

\Mabcd\ 



1 + 5 a b 1 + 5 cd 



(2.3) 



where |-M a f> c d| 2 is the colour- averaged matrix element squared and f a ^ are the usual 
parton distribution functions for flavours a and b. Notice that the coupling constant 
has been taken out of the matrix element for future convenience. Since we are 
interested in the gap cross-section, it is useful to express the kinematics in terms of 

V = Vc + Vd , 

Ay = y c - y d (2.4) 

and define 



lo = x a /x b = e y 

z = x a Xb = pcosh 2 . (2.5) 



2 Ay 
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We choose to define the kinematics for the generic process, 



a(p a , i) b(p b , j) -> c(p c , k) d(p d , I) , (2.6) 

such that parton c is always left-moving, while d is always right-moving; i, j, k, I are 
colour indices. Hence we have that y c > yd, which means Ay > 0, or, equivalently: 

0<-t<§, 

^<-«<s, (2.7) 
where the Mandelstam variables are defined by 

s = (Pa+Pb) 2 , 

2 s (, A v 



t = (pa-Pc) = y 1 - tanh 

« = (Pa-^) 2 = -^fl + tanh^) . (2.8) 



At the Born level, the missing half of the phase space is obtained by exchanging 
t and u. However, when resummation effects are taken into account, this issue is 
more complicated and one must consider different anomalous dimension matrices, as 
explained in Section 3. 

In order to obtain the gaps-between-jets cross-section, we need to integrate 
Eq. (2.3) over the total rapidity of the parton system, y. Fortunately the squared 
matrix element does not depend upon y. This is trivial at LO and it remains true 
after the inclusion of the resummation effects. The gap is defined in terms of the 
distance in rapidity between the final-state partons, Ay, and the jet radius R: 

Y = Ay-2R. (2.9) 

The double differential cross-section is thus 

^HA^b = ^(Ay,Q 2 )|A^ cd | 2 -^-^- (2.10) 

dAy dQ z AS ' 1 + d ab 1 + d cd 

a,b,c,a 

and we have introduced the differential parton luminosity: 

£ ab (Ay,Q 2 ) = ±- f dy UJ-zey' 2 ,Q 2 )f b {^e-y/ 2 ,Q 2 ) , (2.11) 

lZ J-y+ 

where the value of z is determined by Eq. (2.5). The integration limits are 

±t = ±rmnUn- z ,y cv A , (2.12) 
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where the value y cut = 2rj — Y — 2R is obtained by requiring that both jets are 
within the calorimeter acceptance rj. Since for both ATLAS and CMS the edge of 
the hadronic calorimeter is at 4.9 units of rapidity, a sensible choice is 

77 = 4.5 for R = 0A. (2.13) 

We also choose Qo = 20 GeV as the veto scale. This is a fairly conservative choice 
in order to avoid contamination from the underlying event. It could be lowered once 
those effects are better understood. 

The relevant flavour combinations of the parton distribution functions are 



Q (1 \z,Q 2 ) 



1 

2~z 



Q (2 \z,Q 2 ) 



Q^(z,Q 2 



Q (4 \z,Q 2 ) 



S(z,Q 2 ) = 



+ q 
1 

2z 



" dy [Qa(V~^ /2 , Q 2 )qa{V~ze-y'\ Q 2 ) 

V + a 

^,Q 2 )q a (V-ze-^,Q 2 )} , 



y 



dyJ2 [qa{^/ 2 ,Q 2 )q b {^ze-y' 2 ,Q 2 ) 

a^b 

+ q a ( v rzey/ 2 ,Q 2 )q b (V~ze-y/ 2 ,Q 2 )} , 

kz f-+ dv S ^^ e ' /2 ' Q 2 )^*-" 2 , Q 2 ) 

~y a 
+ q a (V~zey/ 2 ,Q 2 )q a (^e-y/ 2 ,Q 2 )] , 

1 fV ^^[g a (^/ 2 ,g 2 )g 6 ( v ^e-^ 2 ,g 2 ) 

a^b 



2z 



+ q a (^ey/ 2 ,Q 2 )q b (^e-y/ 2 ,Q 2 )} , 



1 

2~z 



G(z,Q 2 



a 

{q a (Vz~ey /2 , Q 2 ) + q a (^ 2 , Q 2 )) 9 (y^e^ 2 , Q 2 ) 
{q a {V~ze-y' 2 , Q 2 ) + q a (Vz~e-y/ 2 , Q 2 )) g(^ 2 , Q 2 ) 



(2.14) 



The hadronic cross-section is then computed by collecting together the various par- 
tonic contributions, multiplied by the appropriate parton luminosity: 



dV*) 
dy dQ 2 



pna 2 



\M mq \ 2 Q {1) + \M qq/qq/ \ 2 Q 



4S 

\-Mqqqq 



2 0(2) 



+ \M qq ,qq,\ 2 Q^ + \M 
+ (nf\Mggqq\ 2 + \M 



+ {Uf ~ l)\M q qq,q,\ 2 + \Mqq gg \ 

2 S 



2 ) C< S > 



qgqgi 



9999 1 



(2.15) 
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where we have collected together the contributions coming from the subprocesses 
with four quarks and four antiquarks. Henceforth we write the cross-section with 
a superscript (i), which denotes the number of gluons outside the gap. The tilde 
reminds us that we are generally considering resummed squared matrix elements. It 
is the aim of Section 3 to compute such matrix elements in the case of no out-of-gap 
gluons, i.e. i — 0. We consider the case of one gluon outside the gap in Section 4. In 
the following we also omit the double derivative notation, but it is understood that 
we always consider the double differential cross-section. 



3. Zero gluons outside the gap 

The resummation of the global logarithms is achieved by considering the original 
four-parton matrix element, dressed by virtual gluons "in the gap", with transverse 
momenta above Qq. No out-of-gap gluons are included. The formalism is well known 
and it was first applied to gaps-between-jets in [14]. Here we recap the result and 
also address, in some detail, the issue of the t — u symmetrisation, because we have 
found it to be less than clearly explained in the existing literature. We then study 
the effects of resummation at the LHC and compare our results to Herwig+ + . 

The colour basis independent formalism was introduced in [41] and more re- 
cently applied to the gaps-between-jets cross-section in [26]. In this framework the 
resummation of global logarithms is achieved by considering 

\M\ 2 = Um \e-^e^ r \m ), (3.1) 

» c 

where 

i = m^Q) = -\ -r^«.(*r) (3-2) 

and V c is an averaging factor for initial state colour. The vector \m ) represents the 
Born amplitude and the operator T is the soft anomalous dimension: 

T = l -Yt 2 t + mt a ■ t b + \p&{Y, \^y\){t 2 c + t 2 d ) , (3.3) 
where tj is the colour charge of parton i and 

The operator tf describes the mapping from an m dimensional vector space onto an 
m + 1 dimensional one, as a consequence of emitting a soft gluon with colour a: 

|m + l) = ^y;^ft?|m). (3.5) 
. Pi ■ k 
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The operator t\ represents the colour exchanged in the t-channel, i.e. t t — t a + t c = 
t b + t d . Some care is needed however since the flavour of partons c and d changes 
under the t u symmetrization. If we use the notation 12 — > 34 to label the parton 
flavours then 

t 2 = (tx + t 3 ) 2 = (t 2 + t 4 ) 2 (3.6) 
in the case that parton 3 is identified with parton c, and 

t^^ + t^H^ + ta) 2 (3.7) 

in the case that parton 4 is identified with parton c. 

In order to compute cross-sections, one must fix a basis in the vector space 
spanned by |m) and consequently a representation for the colour operators has to 
be computed. We choose to represent the 4-parton amplitude by M.ij±,h where the 
i,j (k,l) are the colour indices of the incoming (outgoing) partons. In the case of 
the radiating parton being an outgoing quark or an incoming antiquark, the colour 
operators are represented by t%, the generators in the fundamental representation. 
In the case of an incoming gluon, they are T£ c = —if£ c and the sign reverses in 
the case of an incoming quark or an outgoing antiquark or gluon. We always work 
in an orthonormal basis, expecting symmetric results for the anomalous dimension 
matrices [42,43]. The expression in Eq. (3.1) then becomes 

\M\ 2 = ^tr \Mte-^e-^M ] = tr (HS) , (3.8) 

where the hard and soft matrices, H = Aio-M-o/V c and S = e~^ rt e~^ r have been 
introduced. In the following sub-sections we analyse in detail the various channels 
that contribute to the cross-section. For each sub-process we fix a basis and compute 
the hard scattering matrix at the Born level (H) and the soft anomalous dimension 
T. We usually fix Ca = N c = 3 and Cf = ^fN^ = |- 

3.1 Anomalous dimension matrices 

We start by considering the scattering of non-identical quarks, 



qq -> qq', 



in the basis: 



c 2 = ^= \ SuSjk - T^Sikdji ) . (3.9) 
The hard scattering matrix is given by 
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The anomalous dimension matrix computed by projecting Eq. (3.3) onto the basis 
specified by Eq. (3.9) is 



13 



2 

3n et 3 

in\/2 2 i 3v 



3 3 



Pjet + %Y 



(3.11; 



where the subscript 13 indicates that we have considered t 2 = (t x + 1 3 ) 2 (in which 
case the basis is a t-channel singlet-octet basis). 

The anomalous dimension which describes the other half of the phase-space is 
given by 

(2p. ct _|_ ly iny/2 _ V2y \ 
iL§_ hv In 3 + IV- is. ' t 3 ' 12 ^ 
3 3 1 3« ct "r" 6 1 3 / 

The basis specified in Eq. (3.9) is then a w-channel basis. The same matrices describe 
the subprocess with four antiquarks. The resummed squared matrix element for non- 
identical quarks (or antiquarks) is then 

\M qq , qq ,\ = ti{H(t,u)S 13 + H(u,t)S u ] . (3.13) 

In the case of identical quark scattering 

qq -> qq 

the hard scattering matrix is 

hm =a41 * (3 - 14) 



with 



Xi 



t 2 + s 2 



X2 = 3 



9 ' 
IT 

s 2 t 2 + s 2 



ut V? 1 ' 

u 2 + s 2 U 2 + s 2 2s 2 

* 3 = ^2- + 9— "3^- (3 ' 15) 

The resummed squared matrix element for the case of four identical quarks (or an- 
tiquarks) is 

\M qqqq \ = hi [H(t,u)S 13 + H(u,t)S u ] = tr [H(t,u)S 13 ] , (3.16) 

as expected for identical final state particles. 

We now consider the scattering processes involving two quarks and two anti- 
quarks, 

qq -> qq, 
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for which we choose the basis: 



1 ( M 1 



c 2 = ^= [SijSki - -dikdji ) . (3.17) 



The hard scattering matrix is 



where 



H(t,u) = -[ 2 ^ 9A ^|, (3.18) 

9 



9 1^X2 X3 



Xi 



t 2 +M 2 
S 2 



It 2 t 2 + It 2 

u 2 + s 2 lt 2 + u 2 2u 2 , oin . 
* 3 = — + 9— "3rf- (3 ' 19) 
The two anomalous dimension matrices in the basis given by Eq. (3.17) are 

(3.20) 



13 



2 -iivV2 
3H et 3 
— t7r\/2 2 I 3 \z" 7i7r 



3 3 



Pjct + |^ 



and 

2 _ j_ 4y -wry/g , \/2 



14 



_ I 4y -my 2 , V2va 

3 Pjct + 3^ — H^-r 

27r\/2 | y2 y 2 _i_ V 7ct 



(3.21) 



3 1 3^ 3^j ct + ^ 6 

The expression for the resumed matrix element is obtained by adding the two con- 
tributions: 

\M qm \ 2 = tr [H{t, u)S 13 + H(u, t)S u ] ■ (3.22) 

For the unequal flavour process qq' —* qq', one must drop the s-channel terms, so 
that what is left is 

Xi = 0, 

X2 = 0, 

Xs = ! ^ 2 L ^ • ( 3 - 23 ) 
For — > g'(f one drops the t-channel contributions, obtaining 

„ t 2 + it 2 

Xl = 2~, 

s z 

" t2 + u 2 

.. It 2 + u 2 , 
XS = ■ (3.24) 
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We now turn our attention to quark-gluon scattering: 

qg -> qg and qg -> qg . 



A suitable basis is 



1 



ci 

c 2 
c 3 



24 
3 




u jlc L ki 5 



1 



(3.25) 



where now are adjoint indices. In such a basis the hard scattering matrix is 

/ 4xi 



H{t,u) 



1 

24 



^fa lop 4v ^ Xs 
V4v/2X2 4^X2 12X3 J 



(3.26) 



where 



Xi 

X2 
X3 



2 - 



1 - 



su ' 

1 1 2 



u 



3-4- 



2 SM st 

1 t 2 



t 2 2su 

The anomalous dimension matrices turn out to be 





(3.27) 



/ 13 n- -I- 3in 



13 



V2i 



\ 



12 net i 4 



2 



12 ^Jet i 2 ' 



(3.28) 



and 



/ !:!/; j( „ • ¥Y 



14 



12' 



V 



3i7r I 13 
4 6 



\/2i-k \/2 y \/5^7r 

2 Z 







\/2i 



13. , 17y V^iTT 

12n°t "T 12 1 4 4 



2 - 1 x 

\/5y 



(3.29) 



2 2 " 4 4 J 12net ^ 12 - 1 / 

The hard scattering matrix and the anomalous dimensions for antiquark-gluon scat- 
tering have the same expressions as the ones stated above if one chooses a new basis, 
obtained from Eq. (3.25) by exchanging i <-> k and j <-» /. The resummed result is 



\M 



qgqg \ 



tr[H(t,u)S 13 + H(u,t)S u ] ■ 



(3.30) 



For the process 



99 -> 99 
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it is convenient to define an s-channel basis 

1 



Cl 

c 2 
c 3 



24 
3 

20 
1 




U<ljc<>kl J 



The hard scattering matrix is then given by 



3 3 
2^X1 10X1 
3 3 



V 2^X2 2V5 X 2 6X3 7 
with 14 = (iV 2 — l) 2 = 64 and the functions Xi are 



t u 

Xi = ~ + t , 
■u t 



X2 = -2 



u - t t 2 -u 2 



X3 = " + . 

The anomalous dimensions are now 



tu 

t u t 2 + u 2 



/ 2 n . i 13 V 

' 3^J ct + T r 



r 



13 



3t7r 
2 



V 



o 

2 1 



2 - 1 

V5y 



r i4 = 





V 



13 v - 3-i7r 
2 



^2 





2 n . i llV 

3 Pjct T 12 i 4 

_Vlv 2 j 

4 2 3^J et ^ 12 2 / 

\ 

In- 4- ^>Y 

3 "jet ^12 4 



V5 



^ 3^jet+l2^/ 



The resummed result is 

I At 

The process 



ggqq I 



tr [H(t,u)S 13 + H(u,t)S 14 \ 



qq -> ^ 

is related to the previous one by time reversal. In the s-channel basis 

1 



Cl 



c 2 



24 
3 

20 
1 




u lkc L ji 1 



the hard scattering matrix is the same as Eq. (3.32) but with V c = N% = 9. As far as 
the anomalous dimension is concerned the only difference with respect to Eq. (3.34) 
is in the coefficient in front of the jet function pj et : 



/ 3 , 13 V 

1 2^jet + ~e Y 



13 



Y 



3in 
2 





_\/5y 3 



2 1 



2 pjct + 



(3.38) 



Similarly 



ii 



V 



iP'i ct ' 6 


V2 



Y 



Then the resummed result is 

f 



'Sin 
2 







2 

4- IZV V$Y 

y|y 3 , 17 



*7 



l-H*?/ = 2 tr n ) 5 i3 + # («, OS'm] = tr [if (t, u)S 13j 



(3.39) 



(3.40) 



as expected for identical final-state partons. 

Finally we need to compute the four-gluon scattering process 



99 ~> 99 



We choose the following basis: 

f 



Cl = 






1 3 


c 2 = 


2^2 5 ,,c 




1 X f f 


c 3 = 


2V2 3 hkcfjlc 



c 5 



2\/5 
1 

3^ 



2 (^-4; - ^a^fe) 



g fikcfjlc 



The hard scattering matrix is 



H(t,u) 



1 

16 



^9X1 9V2 X i 9V2 X 2 -9^3xi \ 

9^X1 18xi 18X2 -9^X1 

8x3 -9V&X2 





9V2 X 2 18X2 




V -9V3xi -9V6X1 -9^X2 27 X i J 



(3.41) 



(3.42) 
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where 







tu 


St 


Xi = 


1 - 


^ ~ 










It 




St 


ut 




X2 = 


u 2 




+^ 




27 


-( 


su 


X3 = 


T 





SU 

t-u ' 

1 1 st 

4 s 2 4 m 2 



We find for the anomalous dimension matrix: 

/ 



9 /V 

2 V st + ut 



1 t^ 

2 su 



(3.43) 



"Pjetl5x5 + 



V 



3j7T 





3%/2»7r 








4 


4 





3y 

2" 1 


3J7T 


3Vl6l7r 





4 


10 




3i7T 









4 


4 


4 












\/l5i7r 


10 


5 












Ay 5in 
^ J 4 


4 


5 



(3.44) 



/ 



We have already shown that in case of identical final state particles there is no need 
to compute T^. The resummed squared matrix element is thus 



\M 



9<J99\ 



ti[H{t,u)S 13 ] . 



(3.45) 



3.2 Phenomenology and comparison to HERWIGH — h 



In this section we compute the gaps-between-cross-section at the LHC, comparing the 
result obtained from the resummation of global logarithms to the one generated by a 
leading order parton shower Monte Carlo event generator, Herwig+-|- 2.3.0 [45-48]. 
In both cases we use the the latest MSTW LO parton distribution functions [44]. 

As a cross-check, we first compare the result of the Born level calculation. The 
results are presented in Fig. 2. The double differential Born cross-section is computed 
using Eq. (2.15) without any resummation and it is compared to that obtained using 
Herwig+4" at the parton level but without any parton showering. The agreement is 
good, although there is a small (10%) difference that arises because the computation 
of the cross-section from Eq. (2.15) uses the value of the strong coupling constant 
provided by the parton distribution function fit, whilst HERWIG++ uses its own 
one-loop coupling. The leading order MSTW fit produces a larger value of a s and 
this explains the difference. 

We next present results after parton showering and jet finding. The jets are 
defined using the SlSCONE [49] algorithm, with a cone radius R = 0.4 and an overlap 
parameter, / = 0.5. The hardest two jets, which define the extremities of the gap, 
are required to have \yi\ < 4.5 and Q > 20 GeV. We impose a veto such that there 
should be no third jet in-between the pair of hardest jets with kx above 20 GeV. 
To assess the impact of hadronisation corrections, Fig. 3 shows the gap cross-section 
as a function of Q, for Y — 3, and as a function of Y, for Q = 500 GeV, at both 
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Figure 2: The double differential Born cross-section (nb) computed using Eq. (2.15) 
without any resummation (solid lines), compared to the one obtained using Herwig++ 
at the parton level (histograms) and without any parton showering. On the left the cross- 
section is plotted as a function of Q, for Y = 3 (upper lines) and Y = 5 (lower lines). On 
the right it is plotted as a function of Y, for Q = 100 GeV (upper lines) and Q = 500 GeV 
(lower lines). 
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Figure 3: The gap cross-section as a function of Q, for Y = 3, and as a function of Y, for 
Q = 500 GeV, at both parton shower level (solid) and after hadronisation (dotted). 

parton shower level and after hadronisation. The effects are modest and henceforth 
we always compare to the Monte Carlo results after parton showering but before 
hadronisation. We remark that corrections from the underlying event (which we 
have not included) ought also to be small due to our conservative choice of veto 
scale. 

In Fig. 4 we show the gap cross-section, normalised to the Born cross-section (i.e. 
the gap fraction), as a function of Q at two different values of Y and in Fig. 5 it is 
presented as a function of Y at two different values of Q. The solid lines represent the 
results of the (zero out-of-gap gluons) resummation; the dashed lines are obtained 
by omitting the iir terms in the soft anomalous dimension matrices. Neglecting these 
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Figure 4: Plot of the gap fraction as a function of Q, for Y = 3 (left) and for Y = 5 
(right). The dashed lines are obtained by dropping the imaginary part of the evolution 
matrices. 




Figure 5: Plot of the gap fraction as a function of Y, for Q = 100 GeV (left) and for 
Q = 500 GeV (right). The dashed lines are obtained by dropping the imaginary part of 
the evolution matrices. 



contributions corresponds to neglecting the contributions from Coulomb gluons. As 
a consequence, the gap fraction is reduced by 7% at Q = 100 GeV and Y = 3 and 
by as much as 50% at Q — 500 GeV and Y — 5. Large corrections from this source 
herald the breakdown of the parton shower approach. This is as expected, for it 
is these terms that are ultimately reponsible for the colour mixing that generates 
colour singlet exchange, which dominates for large enough Y and small enough Q . 

In Fig. 6 we compare the gap cross-section obtained after resummation (the 
line types are as in Fig. 4 and Fig. 5) to that obtained using Herwig++ after 
parton showering (Q is taken to be the mean of the two leading jets). The broad 
agreement is encouraging and indicates that effects such as energy conservation are 
not too disruptive to the resummed calculation. Nevertheless, the histogram ought to 
be compared to the dotted curve rather than the solid one, because Herwig+4- does 
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Figure 6: Plot of the gap cross-section at Y = 3. The resummed results (line types as in 
Fig. 4 and Fig. 5) are compared to the result of Herwig++. 



not include the Coulomb gluon contributions. We do expect the imposition of energy 
conservation to enhance the resummed gap cross-section. In addition, the influence 
of non-global logarithms (some of which are included in the Monte Carlo) should 
act in the opposite direction and push down the dotted curve. Finally, there is the 
fact that the Monte Carlo operates only at leading N c to consider. As anticipated, 
there are very significant corrections to the parton shower Monte Carlos at the larger 
values of Q/Qo- Of course the resummation would benefit from matching to a NLO 
calculation and this should be done before comparing to data. 



4. Super-leading logarithmic contributions 

In this section we shall evaluate the effect of including the super-leading logarithms 
(SLL). In [27], the coefficients of the SLL have been computed order-by-order in 
perturbation theory up to 0(a^), with respect to the Born cross-section. 1 At this 
order in the perturbative expansion one must consider the contributions coming from 
one and two gluons outside the rapidity gap. We start by computing the impact 
of the super-leading contributions on the hadronic cross-section order-by-order in 



1 Morc precisely, the calculation is performed in the high energy limit, considering only the 
i-channcl gluon exchange contributions to the Born cross-section. 



perturbation theory. To this end we define three different i^-factors: 
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(4.1) 

The SLL contributions to the cross-section are the ones computed in [27]: 



5 



r (2) _ i ^born f a A 5 C 2 , v A Q 



where Cj and V>i are colour factors. Notice that we have chosen the notation in 
order to emphasize the factorial denominators coming from the strongly-ordered fcr- 
integrals. We plot these .fT-factors in Fig. 7 as a function of Q, for Y = 3 (left) and 
for Y = 5 (right). The curves on the left (right) are obtained by considering one 
out-of-gap gluon at 0(a A s ) (dotted) and 0{a b s ) (dashed). The two-gluons outside the 
gap contribution is included in the dash-dotted curves. We also plot these ^-factors 
for two values of Q as a function of the rapidity separation Y in Fig. 8. The plot 
on the left is for Q = 100 GeV, while the one on the right is for Q = 500 GeV. 
The different line styles are as in Fig. 7. From these plots it is clear that, for jets 
with rapidity separation Y > 3 and transverse momentum bigger than 200 GeV, the 
inclusion of fixed order SLL contributions leads to sizeable but unstable i^-factors. At 
every order in perturbation theory, for a given number of out-of-gap gluons, the SLL 
contribution acquires an extra power of Y and L; the colour factor also increases. For 
large enough rapidities and transverse momenta, higher order contributions become 
important and their resummation is necessary. Currently the resummation of SLL 
contributions has been carried out only for the case of one gluon outside the gap, 
and we will discuss it in the next section. Figures 7 and 8 offer some hope that 
the impact of the two-or-more gluons outside the gap contribution may be modest, 
since the difference between the i^-factors which include the two-out-of-gap gluon 
contributions (dash-dotted curves) and the ones that do not (dashed curves) is not 
large. 
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Figure 7: Plots of the i^-factors as a function of Q, for Y = 3 on the left and for Y = 5 
on the right. The curves are obtained by considering the one out-of-gap gluon cross-section 
at 0(aj.) (dotted) and 0{a s ) (dashed). The two-gluons outside of the gap contribution is 
included in the dash-dotted curves. The solid line corresponds to the resummation of the 
one out-of-gap gluon contributions and is explained in Section 4.1. 




Y 



Figure 8: Plots of the .fT-factors as a function of Y, for Q = 100 GeV on the left and 
for Q = 500 GeV on the right. The curves are obtained by considering the one out-of-gap 
gluon cross-section at O(af) (dotted) and 0{a^ s ) (dashed). The two-gluons outside the 
gap contribution is included in the dash-dotted curves. The solid line corresponds to the 
resummation of the one out-of-gap gluon contributions and is explained in Section 4.1. 

4.1 Resummation of one gluon outside the gap contribution 

The aim of this section is to resum the SLL contributions that arise as a result of 
allowing one soft gluon outside the rapidity gap. The general framework is described 
in [25] and [26]. One must now consider both real and virtual corrections to the four- 
parton scattering, each dressed with any number of soft (either Coulomb or eikonal) 
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gluons: 

o rQ ji, r 

\Mx\ 2 = — ^oc s {k T ) (n fl + n v ). (4.5) 

7T JQo K T Jout 

Super-leading logarithms originate when the out-of-gap gluon becomes collinear with 
one of the initial state partons. 

Since our aim is the resummation of the SLL, we can work in the limit where 
the gluon is collinear to one of the initial state partons. We also consider only 
leading terms in the rapidity separation Y, both in the hard scattering and anomalous 
dimension matrices. As a result, the expressions for real and virtual emission simplify 
considerably. 

First we extract the logarithm that arises after integrating over the rapidity of 
the out-of-gap gluon, i.e. 

I-M^T = -- [ Q d ^<k T ) fin f) (Qf + O . (4.6) 

In the case that the out-of-gap gluon is virtual, we need the operators to insert 
the out-of-gap gluon (7) and to evolve the four-parton matrix element (i.e. T in 
Eq. (3.3)). In the collinear (and large Y) limit they are 

7 = z = l,2, (4.7) 

r = Wt 2 t + iTrti • t 2 . (4.8) 

The diagonal nature of 7 ensures that the virtual contribution turns out to be the 
same as the four-parton case multiplied by a Casimir: 

^ Oll = ^(^o|t^ rt e-« r |m ). (4.9) 

v c 

This result coincides with our expectation that a soft and collinear emission leads 
to an extra logarithm whose coefficient is the Casimir, multiplied by the lower order 
result. If the DGLAP plus-prescription would work, the real emission case should 
give the same contribution but with opposite sign. Interestingly, this does not occur. 
In particular, the operator which describes the emission of a real and collinear gluon, 
D^, and the five-parton evolution matrix, A, are 

D£ = (hi - h)X , « = 1,2, 
A = \YT\ + zttTx • T 2 . (4.10) 

The cross-section depends only on (hi — h) ■ (hi — h) = — 1 and the real contribution 
is thus 

tfon = _ 1 <mo | e -«fc rj O)rt t «t e -«Oo,MAt e -e(Qo,fcr)A t o e -f(* r ,0)r| mo) (4n) 
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The non-commutativity of T 2 and z7i"Ti • T 2 prevents us combining the exponentials 
and means that the contribution to the cross-section from real emissions is no longer 
proportional to a Casimir. At fourth order relative to the Born term, this triggers 
the breakdown of the plus-prescription and generates the super-leading logarithms. 

The five-parton evolution matrices (A) have been computed in [42] for qq — > qqg 
and more recently in [50] for the other cases. We have repeated the calculation but 
do not show the explicit results here because the anomalous dimension matrices can 
be quite big, e.g. in the case of five gluons we have a 22 x 22 matrix [50]. Before 
presenting the results of the resummation, some comments are due regarding the 
operator T 2 . As in the four-parton case, it represents the colour exchanged in the 
t-channel and we use 

T 2 = (T x + T 3 + T fc ) 2 = (T 2 + T 4 ) 2 , (4.12) 

when the gluon, denoted by the index k, is collinear to parton 1, while 

T 2 = (T 2 + T 4 + T fc ) 2 = (T x + T 3 ) 2 , (4.13) 

when the gluon is collinear to parton 2. As previously discussed, in order to imple- 
ment the t «-> u symmetrisation these should be replaced by 

Tf = (Tx + T 4 + T fc ) 2 = (T 2 + T 3 ) 2 , (4.14) 

and 

T 2 = (T 2 + T 3 + T k f = (Tx + T 4 ) 2 . (4.15) 

However, after swapping t <-»• u, the hard scattering matrices are exponentially sup- 
pressed in Y and we can neglect them. Once a matrix representation for the different 
operators has been found, Eq. (4.11) becomes 



Q™ 11 = -tr 



(4.16) 



We checked our results by expanding the exponentials in powers of a s and comparing 
them with the fixed order results of [27] up to 0(a b s ). In [26] it was noticed that the 
coefficient of the SLL at 0(a A s ) was the same regardless of whether the out-of-gap 
gluon is collinear to parton 1 or 2. We find that this is no longer true for higher 
orders if the initial state partons are different. 

Once the real and virtual contributions have been computed, for all sub-processes, 
it is straightforward to obtain the hadronic cross-section. We quantify the effect of 
resummation via a i^-factor, i.e. 

, , <t(0) + 0.(1) 

K*> = (4.17) 

The results are added (as the solid lines) to the plots in Fig. 7 and Fig. 8. Gener- 
ally the effects of the SLL are modest, reaching as much as 15% only for jets with 
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Q > 500 GeV and rapidity separations Y > 5. Remember that we have fixed the 
value of the veto scale Q = 20 GeV and that the impact will be more pronounced 
if the veto scale is lowered. Finally, we ought to explain a particular feature of the 
resummed results in comparison to fixed order. In Fig. 7, the resummed curves lie 
closer to the lower order results than to the O (ccf) results, even at low Q, and in 
Fig. 8, for small values of Y, the resummed curves dip below the O (a A s ) curves. This 
behaviour is completely explained by the running of the coupling. The fixed order 
results are as in Eq. (4.2), with a s = a s (Q), while in the resummed case the running 
coupling is evaluated at a s = a s (fcr), i.e. it sits inside the transverse momentum 
integration. 

5. Conclusions and Outlook 

In this paper we have studied the effects of soft gluon resummation on the imple- 
mentation of a jet veto in the region between the hardest pair of jets in dijet events 
at the LHC The bulk of the leading logarithms result of vetoing primary 

wide angle soft gluon emissions in the inter-jet region and we compare the result of a 
direct resummation of these with the HERWIG++ parton shower Monte Carlo. The 
agreement is generally good, however there are significant contributions arising from 
the exchange of Coulomb gluons, especially at large Q/Qo and/or large Y. Such cor- 
rections are not implemented in the parton shower Monte Carlos. Moreover, there is 
a need to improve the resummed results by matching to the fixed order calculation 
at NLO. These observations will have an impact on jet vetoing in Higgs-plus-two-jet 
studies at the LHC. 

We also attempted to quantify the size of the super-leading logarithms that 
occur because gluon emissions that are collinear to one of the incoming hard partons 
are forbidden from radiating back into the veto region. Although only a partial 
calculation can be performed, we did manage to sum the super-leading logarithms 
to all orders in the case of one collinear emission. We found that the effects are 
generally modest but can reach 15% for values of Q/Qo and Y that are well within 
the reach of the LHC. 

Although our attention has focussed on the role and resummation of wide-angle 
soft gluons, there remains much to be done before one can claim to have properly 
calculated the cross-section for "gaps between jets". In particular, it remains to 
match the resummed results to fixed order calculations at NLO, to systematically 
include the non-global logarithmic corrections and to include the high-energy (BFKL) 
logarithms at large Y (a start in this direction was made in [51]). Moreover, as 
discussed in [52,53], the correct implementation of the jet algorithm even corrects 
the Sudakov exponentiation of the primary emissions, e.g. it has been shown that 
fct-clustering introduces an independent correction and similar effects are expected 
using other algorithms. 
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